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Abstract 



, We derive the linearly perturbed matching conditions between a Schwarzschild space- 

rs) ' time region with stationary and axially symmetric perturbations and a FLRW spacetime 

with arbitrary perturbations. The matching hypersurface is also perturbed arbitrarily 
Q ' and, in all cases, the perturbations are decomposed into scalars using the Hodge operator 

on the sphere. This allows us to write down the matching conditions in a compact way. In 
^ • particular, we find that the existence of a perturbed (rotating, stationary and vacuum) 

bi), Schwarzschild cavity in a perturbed FLRW universe forces the cosmological perturba- 

tions to satisfy constraints that link rotational and gravitational wave perturbations. We 
also prove that if the perturbation on the FLRW side vanishes identically, then the vac- 
^ ' uole must be perturbatively static and hence Schwarzschild. By the dual nature of the 

On I problem, the first result translates into links between rotational and gravitational wave 

perturbations on a perturbed Oppenheimer-Snyder model, where the perturbed FLRW 
dust collapses in a perturbed Schwarzschild environment which rotates in equilibrium. 
The second result implies in particular that no region described by FLRW can be a 
source of the Kerr metric. 



1 Introduction 



" I 

^ . A long standing question in cosmology concerns the way large scale dynamics influences the 

^ I behaviour on smaller scales. The most common view regarding this question is that the 

influence of the cosmic expansion on local physics is zero or negligible. The main argument 
supporting this conclusion is based on the Einstein-Straus model [l] which consists of a 
vacuum spherical cavity (described by the Schwarzschild metric, hence static) embedded in 
an expanding dust Friedmann-Lemaitre-Robertson- Walker (FLRW) model. The matching 
between both spacetimes is performed across a timelike hypersurface using the standard 
matching theory in general relativity, usually known as the Darmois matching conditions (we 
refer to [2] for a full account on matching general hypersurfaces). In this model the local 
physics occurs inside the Schwarzschild vacuole which, being static, perceives no effect of the 
cosmological expansion. 

Despite its clear physical interpretation, this model presents serious problems and involves 
a number of idealisations which include a spatially homogeneous and isotropic cosmological 
model, the assumption of spherical symmetry both for the metric inside the vacuole and 
for its boundary, and the assumption of an exact static vacuum in the interior. Indeed, 
more sophisticated models have been constructed by using Lemaitre-Tolman-Bondi regions 
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for the cosmological part (see [21 H]), and other type of cavities O [6], but all these models 
are spherically symmetric. An important question is whether the Einstein-Straus' conclusion 
is robust with respect to various plausible (non-spherically symmetric) generalisations. 

A general property of the matching theory is that solving a matching problem always 
gives, as an immediate consequence, that the "complementary" matching also holds [6]. 
Consequently, the 'a priori' interior and exterior roles assigned to the Schwarzschild and 
FLRW regions in the Einstein-Straus model can be interchanged, and all the previous results 
also apply to the Oppenheimer-Snyder model of collapse [7], in which a spherically symmetric 
FLRW region of dust collapses in a Schwarzschild exterior geometry. 

The first attempts to provide non-spherically symmetric generalisations can be found 
in the works by Cocke in [8] and Shaver & Lake in [9]. The first conclusive discussion 
of different metrics and shapes of the static region appeared in a paper by Senovilla & 
Vera [10], where it was shown that a locally cylindrically symmetric static region cannot be 
matched to an expanding FLRW model across a non-spacelike hypersurface preserving the 
cylindrical symmetry, irrespective of the matter content in the cylindrically symmetric region. 
Therefore, a detailed analysis became necessary in order to decide whether the assumption 
of spherical symmetry of the static region was a fundamental ingredient for the models. This 
analysis was performed by Mars, who showed in two steps, [11] and [12], that a static region 
matched to a FRLW cosmological model is forced to be spherically symmetric (both in shape 
and spacetime geometry) under very weak conditions on the matter content, which include 
vacuum as particular case. Therefore, the only static vacuum region that can be matched to a 
non-static FLRW is an (either interior or exterior) spherically shaped region of Schwarzschild, 
which thus leads to the Einstein-Straus or the Openheimer-Snyder models. 

To summarise, the Einstein-Straus model, consisting of a static vacuole embedded in an 
exact FLRW geometry (necessarily of dust), does not allow any non-spherical generalisation 
and, in this sense, it is unstable. The same applies to the Oppenheimer-Snyder model. 

Regarding the assumption of staticity, it has recently been shown by Nolan &; Vera [13] 
that if a stationary and axially symmetric region is to be matched to a non-static FLRW 
region across a hypersurface preserving the axial symmetry, then the stationary region must 
be static. Hence, the results in [T2] can be applied to any stationary and axisymmetric 
region. In particular, the only stationary and axisymmetric vacuum region that can be matched 
to FLRW preserving the axial symmetry is an ( either interior or exterior) spherically shaped 
region of Schwarzschild. This implies that the Einstein-Straus and Oppenheimer-Snyder 
models cannot be generalized by including stationary rotation in the non-FLRW region. 
Note that this shows in particular that no FLRW axially symmetric region can be a source 
of Kerr. 

Now, two possibilities in trying to generalise such models can be considered: the first is to 
take non-spherical exact solutions which generalise the FLRW region (such as Bianchi mod- 
els) and the second to consider non-spherical perturbations of FLRW. The first possibility was 
considered by Mena, Tavakol & Vera in [TJ], where they studied a matching preserving the 
symmetry of a cylindrically symmetric interior spacetime with locally rotationally symmetric 
spatially homogeneous (but anisotropic) exteriors. This matching resulted in restrictive gen- 
eralisations of the Einstein-Straus model |14j . none of them physically admissible. Therefore 
it is interesting to consider the second possibility, i.e. to construct a perturbed model. 

Perturbed matching conditions have been applied many times in the past: Hartle [15] 
studied first and second order stationary and axisymmetric rigidly rotating perturbations of 
static perfect-fluid balls in vacuum. Chamorro performed the first order matching of a Kerr 
cavity in an expanding perturbed FLRW model [16]. For spherical symmetry, the linearised 
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matching conditions in an arbitrary gauge were studied by Gerlach & Sengupta [T71 [18] and 
by Martm-Garci'a & Gundlach [T^]. There are also studies by Cunningham, Price & Moncrief 
where axial [20] and polar [21] perturbations of the Oppenheimer-Snyder model of collapse 
were derived. However, a matching perturbation theory in general relativity has only recently 
been developed in its full generality for first order [22\ [23] and second order perturbations 
[24] . A critical review about the study of linear perturbations of matched spacetimes including 
gauge problems has been recently presented in [25] . 

Another interesting approach has been followed by Dolezel, Bicak and Deruelle [26], 
who have studied slowly rotating voids in cosmology with a model consisting on an interior 
Minkowskian void, a matter shell between the void and the cosmological model and a FLRW 
universe with a particular type of perturbation describing rotation. We emphasise that, in 
this paper, we focus on generalisations of the Einstein-Straus and the Openheimer-Snyder 
models without surface layers of matter (i.e. such that the Darmois matching conditions are 
satisfied), and we do not restrict the FLRW perturbations in any way. 

In this paper we consider first order (linear) perturbations of the Einstein-Straus model, 
as well as perturbations of the Oppenheimer-Snyder model (as the perturbative matching 
satisfies the same dual property as the full matching). Our background model consists then 
of a Schwarzschild region matched to a FLRW dust cosmological model across a spherically 
symmetric (timelike) hypersurface. We perturb the Schwarzschild part with vacuum sta- 
tionary axially symmetric perturbations and we perturb the matching surface as well as the 
exterior FLRW region with arbitrary perturbations, not restrained to any material content. 

Our approach to this problem consists in exploiting the underlying spherical symmetry of 
the background and of the matching hypersurface as much as possible. This is normally done 
in the literature by resorting to decompositions of all objects in terms of scalar, vector and 
tensor harmonics on the sphere. Our aim is to use an alternative method based on the Hodge 
decomposition of all tensor objects on the sphere in terms of scalars. The two approaches 
are obviously related to each other. However, by working with Hodge scalars we avoid the 
need to deal with infinite series of objects (one for each / and m in the spherical harmonic 
decomposition). In particular, our set of matching conditions has a finite number of equations 
(involving scalars that depend on the three coordinates in the matching hypersurface) instead 
of an infinite collection of matching conditions for functions of only one variable (the time 
coordinate on the matching hypersurface). The equations are therefore much more compact. 

In fact, even when working with 5^ scalars the length of the equations can grow substan- 
tially depending on how they are combined and written down. This is partly due to the level 
of generality we leave for the matching hypersurface and the FLRW perturbations. We have 
taken the effort to combine the equations and to group several terms in each equation in 
such a way that the set of equations becomes reasonably short and manageable. We empha- 
size that the gauge in the FLRW part will be left completely free, so that the readers may 
choose their favourite one. As an example, we rewrite the equations in the Poisson gauge in 
Appendix [Cl 

Deriving and writing down the linearised matching conditions for our perturbation of the 
Einstein-Straus (and Oppenheimer-Snyder) model is the main result of this paper. A more 
detailed analysis of the resulting set is postponed to a later paper. However, in order to show 
the usefulness and power of the equations, we present two applications. The first one is based 
on the observation that the linearised matching conditions can be combined in such a way that 
two equations involving only terms on the FLRW side hold. These equations are therefore 
constraints on the perturbations in the FLRW part that must be necessarily satisfied if an 
interior stationary and axisymmetric perturbed vacuole is present (and hence the local physics 
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can remain unaffected by the cosmic expansion, as generally believed). These two constraints 
link the vector and tensor (linear) FLRW perturbations on the boundary of the region, and 
they imply, basically, that if the perturbed FLRW region contains vector modes (with I > 2 
harmonics) on the boundary, then it must also contain tensor modes. In other words, if 
a (perturbed) FLRW region contains rotational perturbations that reach the boundary of 
a perturbed stationary and axially symmetric vacuum region, then the cosmological model 
there must also carry gravitational waves. Given the estimates of cosmological rotational 
perturbations through observations (see e.g. [27j and references therein), the constraints due 
to the possible existence of stationary and axisymmetric vacuoles would provide estimates of 
cosmological gravitational waves. 

As a second application, we consider the case when the FLRW part of the spacetime 
remains exact (i.e. all perturbations vanish there). As discussed above, the results in [13] 
combined with [12] imply that the interior has to be exactly Schwarzschild provided the match- 
ing hypersurface is assumed to be axially symmetric. Since we allow for non-axially symmetric 
perturbations of the matching hypersurface we can address the question of whether this result 
generalizes to arbitrary hypersurfaces (at the linear level, of course). Our conclusion is that 
indeed this is the case. This result points, once again, into the fragility of the Einstein-Straus 
model against any reasonable generalisation. From the Oppenheimer-Snyder model point of 
view, this also means that a body modeled by a dynamical FLRW model, irrespective of 
its shape and its relative rotation with the exterior, cannot be the source of any stationary 
(non-static), axially symmetric vacuum exterior, in particular of the Kerr metric. 

The paper is organised as follows. In Section [2] we give a very brief summary of the 
linearised matching theory, where we fix our notation. In Section [3] we summarize the stan- 
dard theory of the Hodge decomposition of vectors and symmetric tensors on the sphere and 
apply the theory to first order perturbations of a spherical background in general. The use of 
Hodge decompositions introduces a so-called kernel freedom which is discussed and analysed. 
Section S] is devoted to obtaining the explicit form of vacuum, stationary and axially sym- 
metric perturbations of the Schwarzschild region. After perturbing the background matching 
hypersurface arbitrarily, we obtain, in subsection 14. 3^ the explicit form of the Hodge scalars 
of the perturbed first and second fundamental forms of the matching hypersurface. The 
same procedure is followed in Section [5] for the FLRW side. Section [6] is devoted to obtaining 
the linearised matching conditions for our problem in terms of S*^ scalars. As in any Hodge 
decomposition, the set of equations decompose into odd and even equations, and both sets 
are carefully combined and rewritten to make them as compact as possible. The equations 
given in this section constitute the main result of this paper, so we summarize the hypotheses 
and conclusion in Theorem 16.11 Sections [7] and [8] contain the two applications we present in 
this paper. Section [7] is devoted to deriving the constraints in the FLRW side and Section [8] 
to the matching with an exact FLRW. The paper contains four Appendices. Appendices [Al 
and [B] contain, respectively for the Schwarzschild and the FLRW parts, the full expressions 
for the linear perturbations of the first and second fundamental forms prior to their Hodge 
decomposition. AppendixO specifies our general matching conditions in the Poisson gauge in 
the FLRW part. Finally, Appendix [D] is devoted to the comparison between the expressions 
used in our formalism and the doubly gauge invariant quantities in |23j and |25j . 

Lower case Latin indices at the beginning of the alphabet a, b, ... = 1, 2, 3 refer to tensors 
on the constant cosmic time hypersurface in FLRW, at the middle of the alphabet, i,j, ... = 
1,2,3 are used for tensors on the matching hypersurface. The first upper case Latin indices 
A, B, ... = 2, 3 denote tensors on the sphere while middle indices I, J, ... = 0, 1 are used for 
tensors in the surfaces orthogonal to the spherical orbits. Finally, Greek indices a,/3, ... = 



4 



0, 1, 2, 3 refer to general spacetime tensors. 

2 Linearised perturbed matching theory in brief 

The linearised matching involves perturbing a background which is already constructed from 
the matching of two regions {Mq , g^^^^) and {Mq , g^^^~), with corresponding boundaries 
which are diffeomorphic to each other (thence identified as the matching hypersurface). Tak- 
ing local coordinates on the matching hypersurface amounts to write down two embeddings 

^± : So Mo± (1) 

such that ^>±(5]o) = S^, where are arbitrary coordinates on each of the background 
regions. The coordinate vectors d^i intrinsic to Sq define tangent vectors to the boundaries 

ef^ = Assuming Sq not to be null at any point, there is a unique up to orientation unit 

normal vector n|^^°. The orientation of one of the normals can be chosen arbitrarily but the 
other must be chosen accordingly so that both point to the same region after the matching. 
The first and second fundamental forms are, respectively, q^^^fj = "c^^^^^^^^l^^i, ^'•^^^ = 

— nf^^e^^V^e^^l and the background matching conditions requir^ 

q% = q^%, ^t = ^~y (2) 

Consider now a perturbation of the background metric g^^t ~ 9^^^^ + 9^^^^ and of the 
boundaries via the vector fields = Q ^n^l^ +f^\^±, where are tangent to S^. 
The linearised matching conditions are derived in [2^ and [23] (see also [24] for second order 
matching), and read 

<Pt,=<Plr ^t,=^ip (3) 
with (for a timelike matching hypersurface) 

= C^.q^% + 2Q^k^%+ef-efg^^^tM, (4) 

+ ^.««%ni°)«ni°)/^fc(0)S. - '^ef-ef I,., (5) 

where Di is the three dimensional covariant derivative of {'E,o,q^^^fj) and 

c(l)±o ^ 1 (l)±a y± (1)±„ _y±a (1)± \ 

The tensors q^^^^ and k^^^^ are spacetime gauge invariant by construction, since they are 
objects intrinsically defined on S^, and therefore conditions ([3]) are spacetime gauge invariant. 
Moreover, it turns out that the equations ^ are also hypersurface gauge invariant provided 
the background is properly matched (i.e. once ([2]) hold). 

^ As mentioned earlier, we are not interested in a resulting spacetime with a non-vanishing energy- 
momentum tensor with support on the matching hypersurface (a "shell"), and therefore we do not admit 
jumps in the second fundamental form. 
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The quantities and are unknown a priori, and fulfilling the matching conditions 
requires showing that two vectors exist such that are satisfied. The spacetime gauge 
freedom can be used to fix either or both vectors Z, but these choices have to be avoided, or 
carefully analysed, if additional spacetime gauge choices are made, as otherwise the possible 
matchings might be restricted artificially. On the other hand, the hypersurface gauge can be 
used in order to fix one of the vectors or T~ , but not both (see [25] for a full discussion). 

3 Hodge decomposition of the linearised perturbed matching 
in spherical symmetry 

From now onwards we shall concentrate on spherically symmetric background configurations. 
In order to write the matching conditions in a way which exploits the symmetries we shall use 
the Hodge decomposition on the sphere. There are some good reasons to do that. As outlined 
in the Introduction, the equations naturally inherit the spherical symmetry of the background 
configuration, and thus one expects that any approach based on spherical decompositions 
will render the equations in a simpler form. Previously in the literature, this has been 
implemented by decomposing the relevant quantities into scalar, vector and tensor harmonics 
(see e.g. [23]). Prom a formal point of view, those decompositions are very useful since 
they provide independent sets of equations for the different spectral values, say I and m. 
However, in practice, the harmonic decomposition may become a problem on its own when 
studying explicit models. Moreover, even if the problem can be formally solved for the infinite 
spectrum, the sum convergence of the resulting decomposition should be eventually ensured. 

Using the Hodge decomposition has the advantage that one works with all the different 
{l,m) harmonics of a given quantity at once. In fact, the Hodge scalars correspond in a 
suitable sense, to the resummation of the previous spectral decomposition. Por instance 
instead of using all the (infinite number of) equations that correspond to each value of I and 
m in the matching equations involving (o"(L)o)zm (see [23]), only one equation for the whole 
sum J2imi(^(L)o)imyr (= ^ ^ere) is needed. 

It is clear that one can always go from the Hodge scalars to the spherical harmonics 
decomposition in a straightforward way. However, it is not always easy to rewrite the infinite 
number of expressions appearing in a spectral decomposition in terms of Hodge scalars. 
Working with Hodge scalars involves a finite number of equations and obviously there arise 
no convergence problems (although then one often has to deal with PDP in 3+1 dimensions 
instead of 1+1 equations, which are simpler). Furthermore, their calculation entails a quite 
straightforward procedure, more easily implemented in algebraic computing. We devote the 
Appendix [D] at the end of the paper to relate the functions of the Hodge decomposition used 
in the present paper and the coefficients used in [231 ES] for the scalar, vector and tensor 
harmonic decomposition. 

3.1 The Hodge decomposition on the sphere 

We recall that the Hodge decomposition on tells us that any one-form V on 5^ can be 
canonically decomposed as 

V = dF + *dG, (6) 

where and F and G are functions on the sphere and {'kdG)A = f]^ aDcG is the Hodge 
dual with respect to the round unit metric hAsdx^dx^ = d'd'^ + s\t? 'ddip^ on S"^ . The 
corresponding volume form and covariant derivative are denoted respectively by r]AB and 
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Da. Latin indices A,B,... are raised and lowered with h^s and the orientation is chosen 
so that r]^^ > 0. Furthermore, any symmetric tensor T^b on the sphere can be canonically 
decomposed as 

Tab =DaUb + DbUa + HHab, 

for some 1-form Ua on S'^, which, in turn, can be decomposed as Ua = DaP + {•kdR)A- 

The Hodge decomposition on the sphere has a non trivial kernel, i.e. the zero vector 
and the zero tensor on S"^ can be decomposed in terms non-vanishing scalars, albeit of a 
very special form. First, we consider the kernel corresponding to the vanishing vector on the 
sphere: 

DaF + ti^aDbG = 0. 

Since the only harmonic functions on the sphere are the constant functions, it follows that F 
and G must be independent of the angular coordinates {a;'^}(= {'d,(p}). Regarding the zero 
symmetric tensor Tab, we must solve 

DaUb + DbUa + HhAB = 0, (7) 

which states that U"^ is a conformal Killing vector on the sphere. There are six conformal 
Killing vectors on the sphere: three proper ones and three Killing vectors. They correspond 
to the usual longitudinal and transverse / = 1 vector harmonics respectively (denoted as 
^(i,)^ and V(r)^ with / = 1 in [231 US). Explicit expressions for the conformal Killing vectors 
are obtained using the / = 1 spherical harmonics {m = 1, 2, 3) 

Yi = cos "d, Yi = sin cos (p, Yi = sin -i? sin (p. 

The gradients DaYi^ correspond to three linearly independent proper conformal Killing 
vectors on S^, and their Hodge duals r]^ aDbYi^ correspond to three linearly independent 
Killing vectors on S"^. Therefore, decomposing further Ua as Ua = DaP + rj^ aDbR we 
obtain 

p = Po+E^-^r 

m 

R = -Ro + E RmX^ 

m 

for some eight free coefficients Pqi Pm-, Ro, Pm independent of {??, tp}. Substitution into dTj) 
leads to 

H = 2j2P^Yr. 

m 

3.2 Perturbed matching conditions in terms of scalars 

Our background configuration is spherically symmetric and composed of two spherically sym- 
metric spacetimes {M^ , g'^^'>^) matched across spherically symmetric timelike boundaries 
diffeomorphic to each other. 

Let us for definiteness concentrate on the (iW^, 5^'^)+) spacetime and drop the -|- subindex 
(analogous expressions obviously hold for the {Mq ,g'^^^~) spacetime region). We choose 
coordinates adapted to the spherical symmetry, so that 

= ojijdx^dx-^ + r^{x^){de^ + sin^ ed(f>^), 
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where w/j is a Lorentzian two-dimensional metric and r(x^) > 0. A general spherically 
symmetric boundary can be described by the embedding (or parametric form) 

So:={xO = c^?o)(A),x' = $(O)(A),0 = ^,0 = V5}, (8) 

where {^*} = {A, t?, cp} is a coordinate system in Sq adapted to the spherical symmetry. 
The coordinate tangent vectors to So read 

ex = ^°(o)^xO + ^o)dxiko^ e^ = de\Eo, = 9(^|so, (9) 

where the dot denotes a derivative w.r.t. A. Defining A^^ = —e\^e\'^u!jj\j^Q the unit normal 
to the boundary reads 

n(°) = ^E^(^_^^^dx' + $o^)dx^)b,. (10) 

The sign of A^ corresponds to the choice of orientation of the normal. The first and second 
fundamental forms on Sq read 

q^^\jdCde = -N^dX^ + r2|s„(dt?^ + sin^ M^"^), (11) 
k^^\jdCdi^ = N^ICdX^ + r^\j:Md^^ + sin^ (12) 

where /C = N-'^ex^ ex^^V inf\j^^X = n^^'^^d^i lnr\^^. 

Applying these expressions to each spacetime region (M^ , g^^^^), the background match- 
ing conditions ([2]) become 

Nl = Nl, r+ls, =r_|so, /C+ = /C_, = (13) 

These equations involve scalars on the sphere and therefore do not require any further 
Hodge decomposition. Notice that the embeddings have been implicitly chosen so that 0+ = 
0^{= -d) and (p^ = (j)~{= ip) on the matching hypersurface. This can in principle be modified 
by an arbitrary rigid rotation, which is an intrinsic freedom of any matching preserving the 
symmetry [28]. At the background level, this rigid rotation is irrelevant and can be re- 
absorbed by a coordinate change. However, its effect is not so trivial at the perturbed level 
(see [25] for a discussion on its consequences). 

Consider now an arbitrary linear perturbation and use the Hodge decomposition applied 
to the perturbed first and second fundamental forms. More specifically, we write 

= DAFl + {*dGl)A 
A;«±^ = DAF^ + {*dGl)A 

= DA{DBPl + {*dRl)B)+DB{DAP^ + {*dRl)A)+HlhAB (14) 

= Da (DbPI + {-^dRi)B) + Db [DaPI + {-^dR^A) + H^hAB 

where F^, G%, P^, R%, H^, F^, P^, R\, H^, are scalar functions on 5^ that depend 
on the parameter A. The linearised matching conditions ([3|) can be rewritten as conditions 
involving q^^^^x ^^^^AA together with the functions above. Recalling the existence of a 
non-trivial kernel for the Hodge decomposition, the equalities in ([3]) turn out to be equivalent 
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to 





N^{PS{X)+Pl{X)Y^"^), 



Pi 
Rl 
Hi 



(15) 



Gl 

pk 

Ht 



iVV;5{A), 

iV2(p-(A)+p-{A)y™), 
A^2(;?«(A)+ij-^(A)y™) 
2iV2p«(A)y™, 



where all the functions with a or m subindex (and a different color and size) depend only 
on A and correspond to the kernel freedom discussed above. They will collectively be named 
as kernel functions in what follows. The explicit factor A^^(A)(= = Nl) in front of these 
functions has been added for convenience, as it simplifies some of the expressions below. 

It may seem that adding these kernel functions is redundant, as they do not affect the 
tensors cP''^fj and k'^^'^fy However, it is precisely the fact that we want to impose the matching 
conditions at the level of S'^ scalars that forces us to include them. From a practical point 
of view, the explicit inclusion of the kernel functions in equations (jlSp allows one to choose 
arbitrarily any particular decomposition at either (it) side. In particular, when studying 
existence problems for the matching of two given configurations (decomposed in terms of S"^ 
scalars in an explicit manner) it is important to keep the kernel functions free, as they may 
serve to fulfill conditions which might otherwise seem to be incompatible. For an explicit 
case where the kernel functions turn out to be relevant, see Section [8] below. 

Summarizing, equations (jl5p are the formal linearly perturbed matching conditions writ- 
ten in terms of S"^ scalars. The next task is to evaluate explicitly all the scalars involved 
in the Hodge decomposition of (jH) and ([5]) in the cases we will be considering, namely the 
matching of a Schwarzschild spacetime with a stationary and axially symmetric vacuum linear 
perturbation and a FLRW spacetime with a general linear perturbation. 

4 Perturbed Schwarzschild region (— ) 

First, we describe the perturbations of the Schwarzschild region (denoted by the (— ) sign) 
and derive the perturbed first and second fundamental forms on Sq . 

4.1 Stationary and axially symmetric perturbations of Schwarzschild 

We start by taking a stationary and axially symmetric vacuum metric which in Weyl- 
Papapetrou coordinates can be written as [29] 
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,(0) 



e^^ {dt + Ad(t)f + e 



-2U 



e^^ (dp'^ + dz'^^ + p2rf02j ^ 



(16) 
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where f/, A and k are functions of p and z. Vacuum linear perturbations can be obtained by 
taking derivatives of (jl6p with respect to a perturbation parameter. The result is 



7(1)- 



2U,^2 , o.-2C/„2fc ^ ^(1)^ ^^p2 ^ ^^2^ _ 2e-2f^C/{l)/,2d0^ (17) 



where the perturbation is obviously written in a specific gauge, which we shall denote by 
Weyl gauge. The functions U^^\ A^^^ in (jl7p depend on p and z and satisfy the perturbed 
vacuum equations, written explicitly below. 

The Schwarzschild background is obtained from (jl6p by setting 



1 



2m 



2k 



r {r — 2m) 



(r — m)^ — cos^ 9 ' 



^ = 0, 



where 



p = rsm6'yl , z={r — m)cos( 



so that the background metric reads 
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(0)- 



1 



2m 



dr + 



dr^ 



1 



2m 



+ r2 ((i^^^sin^^-^-^^ 



In these coordinates, the first order perturbation metric (jl7p becomes 

2m ^ 
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(1)- 



-2r^ sin^ 



r'^de'^ 



The perturbed vacuum equations decompose into a pair of decoupled second order PDE for 
;7(i) and A(^) 



,92^7(1) cos^aC/W S^C/W , ,aC/« 



Sr^ sin0 do dO"^ ' ^ dr 

,52^(1) cos 9^(1) (92^(1) , 9^(1) 
r (r — 2m) — ^t-;^ ^ — - — — 1 ^r^^;^ 4m- 



0, 



dr"^ sm9 89 dO"^ dr 
together with a first order system for k^^^ (which is compatible provided (jlSp holds) 



(18) 
(19) 



dr 
d9 



2m sin ^ 



(r — m)2 — m^ cos^ i 

2m sin 9 
(r — m)2 — m^ cos^ i 



r — m sm ( 



dr 



+ cos ( 



d9 



-r (r — 2m) cos 9 — \- {r — m) sin " 



dr 



do 



4.2 Background matching hypersurface 

The general spherically symmetric embedding ^ for Sq is given explicitly by 

So : {t = to(A),r = ro{\),0 = ^,(P = if}, 
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where to{X) and ro(A) are smooth functions {C^ at least) restricted only to the condition that 



Sg is timelike (this implies an upper bound for 



I drp 
dto 



The coordinate tangent vectors ([9]) to 



Sq read now 



to dt + rodr 



do 



'^0 ' 



and the induced metric on Sq is 



q^^dCde 



-NldX^ + rg(A) + sin^ M^'^ 



with 





-1 















The sign of N^{X) will be left free for the moment. The unit normal (jlOj) to 



(20) 



n 



(0) 



1 



2m\ 



1 todr +1 rodi 



ro / 



2m 



ro 



n 



(0) 



— {-rodt + iodr)\^- 



points outwards from the interior Schwarzschild region (increasing r) whenever to > and 
N- > 0. The extrinsic curvature (1121) relative to this normal reads 



k^dcde = ^ 



-toro + toro + 



m 



ro (ro — 2m) r. 



2m 



— ( 1 I UA 



+io {ro - 2m) + sin^ Mip'^ 



which, after comparison with (fT2]) . gives /C_ and /C_ explicitly as 



1 



Smf-Qto 
ro (ro - 2m) 



^(1 



2m 



''0 ) 



-to (ro - 2m) 



(21) 
(22) 



4.3 First order perturbation of the matching hypersurface 

We now derive the perturbed first and second fundamental forms on Sq in terms of scalar 
quantities. We start by considering a general vector 



Z\X, ^, ip)dt + Z\X, ^, ip)dr + Z'^iX, i9, ^)de + Z\X, 1?, ^)d^ 



(23) 



which describes how the matching hypersurface Sq is deformed to first order. Using @ and 
([5|) the perturbed first and second fundamental forms on Sq can be readily computed. The 
results are shown in Appendix O In order to write down the matching conditions in terms of 
scalars, we need to decompose the vector Z'^d^ + Z^d^p (which is tangent to the spherical 
orbits) according to its Hodge decomposition. Explicitly 



dTf 1 dTf 

di!) sin -d dip 



1 dT{ 



+ 



1 d% 



sin -d dip sin -d dd 
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where 7^ (A, tp) and 7^ (A, (^) are 5^ scalars which are defined up to additive functions 
of A. The radial part of Z~ can be also decomposed in the following intrinsic manner 



So 



Q-nL°^+T e 



1 1 



(24) 



where, again, Q^{X,'&,ip) and T~ {\,'d,ip) are scalars on Sq . 

When studying the Hodge decomposition of the first order perturbation tensor S'^^-'^, we 
found it convenient to define two new scalars Q{r, 9) and 'P(r, 9) by 



^(1) 



sm 



d'^V cos9dV 



89^ 



sm ( 



89 



The function Q is defined up to an additive function of r while the kernel freedom in V 
corresponds to Po{r) + Pi{r) cos 9, for arbitrary Pq and Pi. 

As a side remark, we recall that a useful function in any stationary and vacuum spacetime 
is the twist potential Q, which in the axially symmetric case can be written as dA = pe~'^^*dQ, 
(where the Hodge dual operator * refers to the {p,z} plane). The background value of is 
obviously zero. The first order perturbed twist potential Jl^^^ is, in fact, closely related to the 
function Q above. By appropriately restricting the additive function in ^, it can be checked 
that 0(1) = (l - ^) G,r holds (the freedom left in Q is the addition of any function g{r) 

solving (^1 — g^r = c for an arbitrary constant c). 

With all the above definitions, the Hodge decomposition of the angular components of 
q^^^~i and k^^^^j can be computed. The corresponding scalars, following the notation in 



are: 



Fl 

Gl 

Pi 
Rl 

Hi 
Gl 



Rl 



toaiso 1 



2m 



ro'(Tr+P)|Eo 




2r2C/(i) 



2r, 



, cos 9 dV 



° sin 9 89 



(25) 



So 



m 



iofl 



2m 



ro(ro - 2m) 



2mrQ 




2m) tQ 



2m) toTf + ^g 




NlrQ[3m — vq) 
ro(2m - ro) 
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1 



roro 



Af2 



toro - toro - toN_ 2 *o?'o 



2m 



'0 



ro — 2m 



ro — 2m 



2mro 



ro(ro - 2m). 

N'^rQ^m — ro) 



+ 



(2m — ro)^ ro(2m — ro) 



ro - m d (Q 



+ iorlroNZ^ f A^P - 2?7« 



N / a ^ COS0 a^p St/W V 



So 



where A52 denotes the Laplacian on (5^,/ia_b)- 



5 Perturbed FLRW region (+) 

In this section we describe the perturbations of the FLRW region (denoted by a (+) sign) 
and derive the perturbed first and second fundamental forms on the matching hypersurface 

5.1 First order perturbations of FLRW 

On a background FLRW spacetime there exists a coordinate system {t,x'^} in which the 
metric reads 

5(0)+ = a2(r)(-(iT2 + 7,bdx'^(ix''), 

where -iabdx^dx^ = dR^ + f{R, e){de'^ + siiP OdcfP') with / = sinhi?, R, sin i?, for e = -1, 0, 1 
respectively. The covariant derivative associated to ^ab will be denoted by V^. 

Let us now decompose the first order perturbation tensor g^^^"^ into scalar, vector and 
tensor perturbations [30] 

5«+o = -2aH 

Cab) 



a^^^tb = a\-2^^ab + Xa 



with 
where 

and 



Xab = DabX + 2V(„n) + n, 



Dab = VaVb - -7abV^ 



v'Ya = n,-^ = 0, v^n^fe = 0. (26) 

The vector term Wa can be decomposed further into its irreducible parts: 

Wa = daW + Wa 

with 

V^Wa = 0. (27) 



The evolution and constraint equations for each mode, in any gauge, are given e.g. in [3T] and 
can be written in a closed form after a gauge is specified. In this paper we intend to derive 
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perturbed matching conditions which can be used in any (+) spacetime gauge. A gauge will 
only be specified in Appendix O where the Poisson gauge will be chosen as an example of 
how the perturbed matching equations simplify from their general expressions to a specific 
gauge. 

Next, we introduce 5^ scalars Wi, W2, 3^i, 3^2, Qi, 22, ^i, ^2 and Ti according to 
Wede + W^d^ = dWi+*dW2, 

Yede + Y^d<i) = dyi + i.dy2, 

URgde + UR^d4> = dQi+*dQ2, 

Uab = DA{DBUi + {^dU2)B)+DB{DAUi + {^dU2)A)+nhAB. 

where Hab^ * and Da refer here to the coordinates {^,0}. The trace- free condition on Uab 
gives 

iiRR = -^{As2Ui + rL). 

As before, these scalars are defined up the kernel of the Hodge operator, which in this case 
involves functions of r and r. Concretely, each one of Wi, W2, 3^i, 3^2, Qi, Qi, admits the 
freedom Wi — > Wi + wi{t,R), etc., while Ui (and similarly U2) is defined up to Ui — > 
Ui + ui{t, R) + uim{T, R)Y{^ which implies Ti ^ TC + 2uimiT, R)Yi^. Another interpretation 
of this freedom is that Wi,...,Q2 do not contribute to the I = harmonic sector of the 
perturbations, and that Ui,U2,'H do not contribute to the / = 0, 1 sectors, since one can 
always choose the kernels such that >Vi(;=o) = 0> etc... This is made explicit in the relationship 
between these scalar functions and the harmonic decompositions in [231 I25j . as shown in 
Appendix [Pl 

The constraints ()26p and (j27p in terms of the Hodge scalars read 

j^As2yi + '^YR + YR' = 0, 
j^As2Wi + ^WR + W'R = 0, 

As2Qi-2As2Ui' -'^{As2Ui + n)-2n' = 0, 

2As2Ui + 2Ui +n + fQi' + 2ff'Q^-Bi{T,R) = 0, 
A52^/2 + 2^/2 + /'Q2' + 2//'Q2-i32(r,i?) = 0, 

where the prime denotes derivative with respect to R and Bi and B2 are arbitrary functions. 
The latter functions arise because we need to Hodge decompose the equations V^IiaB = 0, 
from which two extra kernel functions appear. 

The relationship between the Hodge scalars and the Mukohyama variables for perturba- 
tions of spherical backgrounds is summarized in Appendix |Dj 

5.2 Background matching hypersurface 

In the Einstein-Straus and Oppenheimer-Snyder models the matching hypersurface is comov- 
ing with respect to the FLRW fiow. In fact, it is now known that this is necessary for any 
matching of a static, or stationary and axisymmetric, vacuum region to a FLRW spacetime 
[T2| I13j. The matching hypersurface Sq" is therefore of the form 

■.{T = x,R = R^,e = ^,4> = ^}, 
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where Rc is a constant. The tangent vectors are 



-^0 



and the first fundamental form is 

q^^dCdi^ = a| [-d\^ + fl (d^"^ + sin^ 

where fc = f{Rc-,e) and oe = a|so = a(A). Comparing this expression with (fTT|) we have 
r^lso = oe/c and N'^ = a|. The unit normal to pointing towards the direction in which 
R increases reads 

n+ = — Or , n+ = asdi^lso, 
and a simple calculation gives the second fundamental form on Sq' to be 

k^^dCdi^ = as/c/^ {d^^ + sin2 ^^d^^) , 
where = f'{R,£)\R=R^- Comparing this expression with (jl2p we find 

/C+ = 0, £+ = /7(as/c). (28) 



5.3 First order perturbation of the matching hypersurface 

The first order perturbation of T,q is defined by a vector field at points on Sq". Similarly 
to the case of the Schwarzschild region, we decompose as 



T+dr + ^dR + - J-^] do + + a 

^ a ^ \ d-d sin-d dip j ^ I sin^ -d dip sin d-d j 



where T+, Q'^, Tj^^ and 7^"^ depend on {X^-d^ip}. The Hodge decomposition of the angular 
parts of q'^^^tj and k^^^f^ (with explicit expressions given in Appendix iBl) in terms of the 5^ 
scalars introduced above can be found after a straightforward but somewhat long calculation. 
Recalling the notation in (1140 . the result i^ 



n 

Pi 



Rl 



Fl 



4 {%^fc+m 



So 



So 



«l(/cV+^2+3^2)^ , 

So 

+2aE (q+fjc + as/c'T^ 



Y -Wr--x' - 3>i' - Qi - 



(29) 



^At some points we slightly abuse the notation and use dot to denote both derivative with respect to r, 
and derivative with respect to A. On the matching hypersurface they obviously coincide as r = A and R, 0, (j) 
do not depend on A there. 
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HI 



+ 



fc 



as 



^ W2'-3^2 -Q2 + 



fc 



2^ 



So 



2 

/c 



aE (/c/c'r2+ + |3^2 + ^^/2'-^Q2 



3^2 + /c%- 



1 A Q- 

7:X 



So 



So 



fc V6 



as 



This concludes the decomposition of the perturbation. Our next aim is to write down and 
discuss the matching conditions. 



6 Matching conditions 

6.1 Background matching conditions: the Einstein-Straus and Oppenheimer- 
Snyder models 

The results in this subsection are well-known, but we reproduce their derivation for com- 
pleteness. The background matching conditions are obtained simply by particularising the 
equations p3p (which correspond to ([2]) in spherical symmetry) to the Schwarzschild region 
and the FLRW region. The second equation in (jl3p implies 

^•o = /cOe- (30) 

Inserting this into (j20p . and using its derivative along A, i.e. ro = /c^E, the first equation in 
([13]) leads to a quadratic equation for namely 

tl = ^- J" . (31) 

(l - fas) 

From N'^ = as we write = cas with cr = ±1 (recall that we want to keep the orientation 
of the normal arbitrary). The fourth equation in (|13p . together with expressions ()22p and 
(j28p give the following linear equation for to 

io = ^^lM_, (32) 

which inserted into (i3T]) yields 

. 2 2 2mas 
as + eag = 3 , 

Jc 
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after using /^^ = 1 — e/^. This ODE for the scale factor is exactly the Friedmann equation 
for dust (restricted to points on Sq), as expected. This is just a consequence of the Israel 
conditions, which impose the equality of certain components of the energy-momentum tensor 
on both sides of the matching hypersurface. 

The last matching condition, namely the third equation in (jl3p . is automatically fulfilled 
once (|30|) . (|32|) and the dust Friedmann equation hold, as a straightforward calculation shows. 

Summarizing, given the necessary condition that the FLRW background is dust, the 
matching conditions are satisfied if and only if (I30p and (I32p hold. We can now study the 
linearly perturbed matching. 



6.2 Linearised matching 

The linearised matching conditions ^ correspond to equating the expressions for q^^^fj and 
given in Appendices [A] and [Bj However, as discussed in Section [3. 2^ the angular com- 
ponents are much better handled if the underlying spherical symmetry is exploited through 
the Hodge decomposition, which allows us to work exclusively in terms of scalars. Thus, 
the full set of matching conditions is given by (|15p after using (j25p and (|29p . together with 
the non-angular expressions (f54l) , (156]) , (fSSl) , and (f57l) , given in Appendices |A] and [B1 

The Hodge decomposition in terms of scalars involves two types of objects depending 
on their behaviour under reflection. For instance, in the decomposition ([6]) the scalar F 
remains unchanged while G changes sign under a reflection. The former scalar is then named 
even while the latter is named odd. This splitting behaviour occurs in any decomposition in 
terms of Hodge potentials. In particular, our linearised matching conditions must split into 
equations involving only even scalars and equations involving only odd scalars. We denote 
them simply as the even and odd sets of equations. 

The odd set is simpler to handle. It is not difficult to see that the equations can be 
rewritten as the following four relations 

T+ + f~^[U2 + y2-{R'o+R'Ly-n] = Tf, (33) 
W2-G;5-^[Z^2+3^2-K+iJWlso] = -^Qfca^\ (34) 

W2' - 2asGS - ^[^2' + 3^2' - 2as(/?,';+i?,f„yr)|So] 

Q2- [Z^2'-2as(i?5+iJ-„yr)] + 2/-Vc[^2-K+i?^n™)] = -aGa^^aj^U (36) 

The even set of equations requires a much more lengthy and subtle analysis. After carefully 
combining the equations, and defining 6Q = Q"*" — Q~ , it turns out that they can be written 
as the following eight equations 

r+ + /-2 [u, + y, + ^x- ips+p^^vn] = r,-+v, (37) 
T+-[Wi + w- -^[Ui+yi + ^x- {ps+p^y^nl eo]] = /' «s ^ + , (38) 



Qi + Yr + x' - \Ui - 2aE{R5 17")] 



+ 2/-Vc[^i - (P^.+P^nvn] + — = -fcfcd^^, (39) 

as or 
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1 d 

as dX 



as 



So 



= 7 n ^ + K^fcO'^f- - + fcaj]{fcaj:e - 2m)— y 



+ 



as/c e - 2m 
as/c - 2m 

Wi -Wr- 2aj:F« 
So 2dE/c/c 



as/c - 2m 



dX 
(2C/(i) 



9r 



[^^i' + yi' - 2aE(/?,';+i?;5,yr)lso] 



dX 



5Q 

OS 



dr 



as aA 



as 



d (bO 



+ 



2/^ase + as/c 



as/c - 2m 
(as/c - 2m)2 



dX \as 
6m 

as/c 



I So 

/c ase 



as as/c /, 



1 



Qj.2 ' 



Os/c 



2m 



5r as/c 
as/^(2C/«-Ae7'), 



2m 



as/^|-AeP 
or 



D 

as/^ 



A(x" + 2<^)-^{Wr + W') 
2jc as 



fc 



+ -^{As2x' - x'(3 - 4e/i)) - (/,: 



Sp J 
/ 



2(1 



2/c^e)^ + 



1 



+ 



+ 



as as/c 
/^(as/c + 2m - 2as/^ ) 

/c(as/c - 2m) 
/c(as/c + as/c e - 4m) 



(2m - 2as/c^e 



2e 



3m 



Q- 



^ as/c 
as as/c 
1 



Wi + W - 
1 



A«P 



9r 

/c cos 3 



as 



2P 



or J 



Qfcf, 



7A52X + 



X" 



So /c 



COS 







sm I 



dd 



[/{!) 



(40) 



(41) 



(42) 



(43) 



(44) 



This set of twelve equations represent the full set of linearised matching conditions for our 
problem. They are valid for any FLRW gauge and any hypersurface gauge. Moreover, they 
include the twenty kernel functions f^^, g^, i-;^, fijj, P,?,, Rl^, F^, G^, P^';, P[j, Pg, R'^ in order to 
allow for any choice of Hodge decomposition on either (±) side. Depending on the problem, 
these kernel functions may play a role. For instance, if the aim is to determine perturbations 
in FLRW given perturbations in the Schwarzschild region, then the kernel functions can be 
put to zero without loss of generality since, in that problem, one is constructing the exterior 
data and changing the kernel functions does not affect the metric perturbations. However, 
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in a situation when two specific perturbations are given and the problem is to determine 
whether they match at the hnear level, then the kernel functions become relevant and cannot 
be dropped a priori. 

The expressions above are written in such a way that the (+) and the (— ) objects are 
kept on the left and right hand sides of the equations, respectively. The only exception being 
the difference 5Q which we found convenient to use with a pivotal role. 

The three equations (f33l) , (f37l) , (pSl) determine the difference vector T'^ — T~ . Recall that 
this difference vector is tangent to the background matching hypersurface and corresponds 
to the freedom in perturbing points within the hypersurface without deforming it as a set of 
points [25j . Recall also that by choosing the appropriate hypersurface gauge one can fix either 
T"*^ or T~ arbitrarily (but not both), and that the difference T"*" — T~ is independent of such 
choice. Thus, the three equations ([331) . (f37|) . (f38]) do not provide any essential information 
concerning the metric perturbations at either (it) side or the shape of the perturbed Sq 
(defined by 5Q). We have been careful in rearranging the remaining equations so that the 
difference — T~ does not appear. So, this somewhat superfluous information gets, in this 
way, separated from the remaining (more relevant) restrictions. 

We summarize the results of this section in the form of a theorem 

Theorem 6.1 Let an Einstein- Straus or Oppenheimer-Snyder spacetime geometry be linearly 
perturbed in such a way that the perturbations inside the Schwarzschild region are stationary, 
axially symmetric and vacuum, while the perturbations of the matching hypersurface and 
of the FLRW region are arbitrary. Assume also that the Weyl gauge has been chosen for 
the Schwarzschild perturbation and that the Hodge decomposition has been used to write all 
tensors on the sphere in terms of scalars. 

Then, the linearised matching conditions are satisfied (and hence a perturbed model is 
obtained) if and only if the equations (3^ , < fg5|) . for the odd part, and the equations 
(3^, ggj, gip, (O), glp, ggi for the even part are fulfilled. 

When a specific gauge is used on the FLRW side, the equations above simplify (sometimes 
notably). As an example, we present in AppendixIClthe linearised matching conditions for the 
particular case of a flat e = FLRW region in the Poisson gauge, for which W = x = ^a = ^■ 

It is also worth noticing that the Einstein equations have been used at the background 
level, but the linearised equations for dust at the FLRW region have not been used anywhere. 
Thus, the equations apply to any perturbation of FLRW regardless of the matter content being 
described. The same comment applies to the Schwarzschild side except for the fact that the 
perturbations have been restricted a priori to being stationary and axially symmetric and 
that the form (fT7|) uses part of the vacuum field equations (in particular, it uses the fact that 
/9 is a flat harmonic function, which allows the metric to be written in the form ()16|) ). The 
reader may have also noticed that while the gauge of the FLRW is kept free, the gauge in the 
Schwarzschild cavity has been fixed from the very beginning. The reason for such a different 
treatment is that we implicitly regard the perturbed Schwarzschild metric as a source for the 
FLRW perturbations, which then become the unknowns. 

Having obtained the equations in a simple and compact form (they may be compared 
with the equations that would result from equating all components in q^^'^^j and k^^^^j in 
Appendix |A] with their pairings q^^^tj ^''■^^ ^^^^tj ™ Appendix [B]) , our aim now is to extract 
some of their direct consequences. A more detailed analysis of the equations is postponed to 
a subsequent paper. 
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7 Constraints on the FLRW side 



The first important general consequence of tfie linearised matcliing equations arises by simply 
considering equations (|34|) and (j36|) . Isolating Q from both equations we arrive at the relation 



= Q2 - W - 2aE(i?s+ij^„ir)] + Vc'f'M - {Ri+Ri^vn], (45) 



which only involves objects on the FLRW side. Therefore, this equation constitutes a con- 
straint on the FLRW perturbations on T,q, irrespective of the (stationary and axisymmetric) 
perturbations of the Schwarzschild region, and links the vector perturbations represented by 
the gauge invariant vector perturbatiorH — dy2/dT and the tensor perturbations driven 
by ^2 and Q2 . 

Note also that although there are kernel terms in the form of 1-7" and Rq+W^y^^, 

which would contribute to the / = 0, 1 harmonics, these can be in principle absorbed into U2 
and U2' respectively and do not affect the value of the tensor perturbations. Thus, equation 
()45p implies that if there are no tensor perturbations, i.e. U2 = Q2 = then, on So, 
W2 — dy2/dT cannot contain harmonics with / > 2. 

In terms of the doubly gauge invariant perturbation variables of Mukohyama [23^ I25j . 
defined for I > 2, the constraint (j45p restricted to ^ > 2 is equivalent to the set of equations 



^/o+ =° -2k+ for all (/ > 2, m), (46) 

as it can be easily checked by using the relations in Appendix iDl together with (I29p . and the 
expressions for the doubly gauge invariants in [23]. 

The main result of this section is then summarised in the following theorem: 

Theorem 7.1 Let a region of a general perturbed dust FLRW be (perturbatively) matched 
across a non-null hypersurface to a region of a stationary and axisymmetric, vacuum per- 
turbation of Schwarzschild. If the perturbed FLRW contains a vector perturbation with I > 2 
harmonics on Sq, then the FLRW region must also contain tensor perturbations on Sq. 

In other words, if the FLRW side contains rotational perturbations, then it must also contain 
gravitational waves, irrespective of the matter content described by the perturbation. 

At points where as 7^ there is a second constraint on the FLRW side which is obtained 
by differentiating along A and using ([35|) to isolate ^,r|Eo = ■^{G\t,o)/ id's fc)- This 
second constraint relates the values on So of W2 — dy2/dT and b(2, with their first and second 
derivatives, and reads 



aE/c{(l-aE/D 



W2'- 1^(^/2' + 3^2' I So) 



W2 - g;; - —[U2 + J2 - (Rl+R'LYnlY^o] (3m + as/c - 2eas/e^) 



-as/^ |W2 - -^{U2 + 3>2|so) - {1^2 + y2')'^^^^—j^-Gl+Rl+Rl.Y-r\ = 0. (47) 



3 



Note that Wa — dYa/dr is the (only) gauge invariant vector linear perturbation [30]. W2 — dy2/dT 
corresponds, then, to the divergence-free and odd part of the vector perturbation in FLRW. 
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8 Matching perturbed Schwarzschild with exact FLRW 



Nolan and Vera [13] have proved that in order to match stationary axially symmetric vacuum 
regions to FLRW regions across hypersurfaces preserving the axial symmetry, the vacuum 
part must be static. As an application of our formalism above, we shall show that their result 
can be generalised to arbitrary matching hypersurfaces (not necessarily axially symmetric) 
to first order in approximation theory. 

Since we want to keep the FLRW exact, we set all the FLRW perturbations equal to zero. 
Notice that this entails a choice of Hodge scalar functions in the FLRW part, and therefore 
the kernel functions in the matching conditions must be kept free. 

Let us start by considering equation ([M]) . which differentiated with respect to and using 
the fact that Q does not depend on (p, leads to R| = i?| = plus 

.en I sin^ 'd ■ , . , 

A'^'^k, = -a^-^Rl. (48) 

Jc 

In order to determine Rl we use the constraints ()45p and (|47p derived in the previous section. 
Setting all FLRW perturbations equal to zero, and extracting the coefficient in of equation 
(ji5]) we get 

Inserting this expression into the coefficient of the second constraint ()47p yields a second 
order ODE for RfiX), namely 

aE(2m - aT.fc)K + aj:{a^fc - 4m)i?« = 0, 

which can be solved to give 

Rl = , , (49) 

where C is an arbitrary integration constant. Thus, ()48p becomes 

^"'i..=-^g ,,;'"f'! V (50) 

/c(«e/c - 2m) 

Our aim is to show that the interior source must be static in this case. Since we are working at 
the perturbative level, first of all we need to determine the necessary and sufficient condition 
that ensures that a given perturbation of a static background remains static to that order 
of approximation. To do that, consider an arbitrary metric ga/3 with a static Killing vector 

^ and a first order perturbation metric g^^ which admits a perturbative Killing vector, i.e. 

there exists a vector such that ^ + is (to first order) a Killing vector of §^(3 + ■ 
Our aim is to find the condition that has to be imposed to ensure that this vector is static 
(to first order). 

Let us, first of all, lower its indices and define |q = (ffa/j + (7^/])('^'^ + ^^^^^) ~ '^"+5a/3^^ + 
(,a^ + 0{2). The perturbed Killing vector is static (to first order) if and only if the linear 
term of C[ad/3^ry] vanishes, or equivalently 

M Ad^ + ^AdM = 0, 
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where M = ^^^^ + g[^\ and g^^^ a = Now, staticity of ^ is equivalent to $ Aw = d$ 

for some one- form w, which allows us to rewrite the linear staticity condition as 

$A{dM-MAw) = 0. (51) 

By construction, this equation gives the necessary and sufficient condition for having a static 
first order perturbation. 

If the static background is moreover axially symmetric and the perturbation is stationary 
and axially symmetric, with no further symmetries, the vector is restricted to having the 
form = + bff, where and a and b are arbitrary constants, and ff is the axial Killing 
vector of the background that remains preserved in the perturbation (i.e. the one fulfilling 

£vM = 0)- 

For a static background of the form (116^ (with A = 0) we have 

^ = Ni:dt, T] = Nrid(f>, 

where and A'^^ are respectively the norms of ^ and fj. For a Schwarzschild background 
they are given by 

iV^ = -(^l-— j, N, = r'sm'6[l-—j . (52) 

Computing d$, = N^^dN^ A ^ we obtain w = —N^^dN^. On the other hand, using (jl7p for 
the metric perturbation, we find 



and thus 



It is now straightforward to rewrite (jSip as 



M=(a + 2[/W)^+(^6 + A«^j rj. 



which simplifies to 



NJ \ iV 



—dNr, - —dN. 



0, 



Integrating and using ()52p we find the expression 

j,3 gjj^2 r 



which restricted onto Sq reads 



A' ' = b h const, 

r — 2m 



Aii)\ j,f!^M^ + const. (53) 
asjc - 2m 

Moreover, the fact that the perturbation A^^^ is time-independent implies that its knowledge 
on the matching hypersurface implies its knowledge on the whole region swept by the 
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range of variation of ro(A). We therefore conclude that any interior perturbation which takes 
the form on the boundary defines a static perturbation, at least on the range of variation 
of ro(A). 

Recalling that both fc and are constants, in view of the expressions (I50p and ()53p we 
have proven: 

Theorem 8.1 The most general stationary and axially symmetric first order vacuum pertur- 
bation of a Schwarzschild metric, matching to an exact FLRW geometry across any linearly 
perturbed non-null surface, must be static on the range of variation ofro{X). 
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A Perturbed (— ) region first and second fundamental forms 



Using the definitions of subsections 14.11 14.21 and 14.31 in expressions ([1|)-(I5]) we find that the 
first order perturbation of the induced metric on Sq has the following components 



.(1)- 



AA 



+ 



'0 



2m 



dX 



to 



ro . dZ^ 



-ro- 



+ 



mZ 



,(1)- 



ro — 2m 

dz^ 



ro-2m'" dX ' {rQ-2mf ° 

(54) 



Ai? 



o(i)- 



y -dip 



dX 



2m\ dZ^ . ro dZ^ 
1 + -T^ro- 



2 ■ 2 ,5^' 



2roZ^ + 2r2 



ro J di) ro — 2m d'd ' 

^ 2m\ (dZ^ \ . 

1 H^+^^'^|eo to + 



ro . dZ^ 



dz^ 



+ 



ro J \dip 



ro- 



ro — 2m dip 



So 



dip 



sin^ -i? 



dZ^ 



,(1)- 



dz^ 



ipip 



2roZ sin + 2ro sini? sin 

t?-^ -C/(^)|so sini9 + Z2 cos 1? 



and that the first order perturbation of the extrinsic curvature has the following components 



AA 



J_ firo-2mf.,dU^ ^ ^^'^-^^^l (k^^^ - W^'^ 
\ r^ or rg V 



2m (ro - 2m) dZ'^ 

n tn 



dX 



-to 



+ ^0 4 -Z^ + tori — + tl^„k^^^ 



9A2 



+ tl 



dr dr 



23 



+ 



iV2 



.2 QzO 



ro — 2m dX tq — 2m dX 



° ° (ro - 2m)2t2 _ rgrS \^ ro (ro - 2m) 



+ 



+ ro- 



mrQ 



( 5m — 6ro -2 „ .2 2ro — 3m 



to + 3r, 



° (ro - 2m) 2 



+ 



iliro - 2m) 



(ro - 2m)2i2 _ rgrg 
2 Srgrg 5^0 



ro (ro - 2m) 



dX 



.^ {ro-2my4+3r^r^ ^^(^^ 
ro (ro - 2m) 



+ 
1 



mtoro 



5ig(ro - 2m) 



2 - Srgrg 



t^(ro - 2m) 2 
-to- 



2^2 



rnr, 



0' 



ro (ro - 2m) 



dX 



So 



1 

+ro 
+ro 



dXd^ 
fd^Z^ 



+ to (?'o - 2m) 



'dz^ 

dX 



mdZ^ .\ . 

-to + toTQ 



dZ^ 



'0 



+ 2to- 



\ (ro - 37n)r§azO^ 



ro — 2m d^d 



So 



\dXd-d do do I ro (ro — 2m) 9?? 

^ .92^1 , o m.^dZ^\ (ro-27n)2.2M« 



(. 



+ 



1 ■ dZ^ (ro — 3m) ro dZ^ \ m (ro — 2m) 



-to 



\dXd(p ro — 2m d(f ro (ro — 2m) d(p 

rodA^ _ (ro - 3m)ro ^^(i)^ 
2 5r ro (ro — 2m) 



So 



1 



d^Z 



5??2 



+ ^0 



ro 



Z^ + to (ro - 2m) 2—- + ro 



dz^ 



dr 



ro- 



dr 



-+ 



+ ro 



(d'Z^\ tororg / dZ' 1 
I 9??2 ) + iv2 1 io 



+ 



2mro 



iV2(ro - 2m) 



^1 _^ ^o'^o^'o 92"^ 



Af2 dX 



So 



1 



+ to- 



cos 



+ro 
1 



5Z0 



+ 



sin-!? d(p 
1 a^(^) COST? / dZ° 



-to 



yd^d'd ' 2 do 
' 92^1 



+ to ('"o - 2m) 



+ sm i)- 



5i9 



sill"!? \ dip 



So 



5y^2 



9^-*^ (ro — m) sin2 i? -, 2mrQ sin2 -dZ^ 
sin 17 cos v—;r7r H 



d^ 



ro 



( dZ^ 
+ (ro - 2m) 2sm2i9^— 



+ 2 cos sin T?Z2 — sin2 i? ro 



N'i{ro - 2m) 

duw 



dr 



ro (d^ 
+ io [ 



+ ^ sin^ + ^ sin2 W A:« - 2[/« 



dX Nl 



cos "!? sin "!? 



1 



9?? 



So 
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B Perturbed (+)— region first and second fundamental forms 

Using the definitions of subsections I5.H 15.21 and 15.31 in expressions ([I])-® we find that the 
first order perturbation of the induced metric on T,q has the following components 



,(1)- 



AA 



y Ai? 



,(!) + 



,(!) + 



9^ ^rrlSo 



2a| 



dT+ 
dT+ 



2aEdsT+, 



(56) 



(1)+| 2 2f2(d%+ 1 d%+\ 

2 ' 



5^'^^|so + 2aE/c/^g+ + 2a|/, 



+ 4/c 



dip 



as 



1 d^T+ ^ COST? 9T2+^ 



9 e4>\^o 



(1) 



1 a2ro^ 



+ 



sin ' sin^ -d dip J 

2 2cosi?5T+' 



sin^ dif"^ sin ?? ' sin d-ddip 

+ 2a^fcfcQ^ sin2 + 2aj,flh^T+ sin^ i9 + 

/'^^'^f , .......^-^i^ . ^^^^^^ 



sin {} dip 



I ; + cost? sin t?^^;^ + sin t?^^-^ „ 
\ a(/9^ ov ovoip dip 



q'^-'2 

cos v- 



and that the first order perturbation of the extrinsic curvature has the following components 
as 9r a| ^ 2as 5i? dX^ 



AA 



Ai9 



as dX 
1 5^^^+. 



+ 



OS _ 

,as 
1 9^'^ 



a2 



(57) 



So 



2as 5i? 2as 86 



Q ^ „2 9 m 



1 dg(^)+ d^Q+ , dsag 



2as (?r 



+ 



dXd^ as 5^9 



^■'"■'^[dXd^ sm^dXdip) 



1 5 



(1)' 



1 9 



1^ .Q.q 



2as 2as 

fcfc 



1 dg(^% d^Q+ , ds^Q 



2os dr 



+ 



dXdip as d(f 



\dXd^ 



2as 



(1)+ «s/c (1)+ , 
y RR ^ „2 9 rR ^ 



So 
1 



2as 5i? 
2 



^dg(^B. 
as 9^ 



Re 



+ 



as d-d^ 



2 2 COS 1? dT^ 
sin ■!? d^dip sin^ i? Sy? 



as dX 



2 

c 



1 



R(t> 



2as \ 9ii 



do 



+ 



COS "!? 

as sini?^' '■^'^ 



So 
(1): 



cos^?aQ^ 



+ 



dddip sin 
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sin^ ^/c/e m+ , 
cos "!? sin "!? 



2— -— ^ + sini? 2 
sin a(/9 

a-s sin^ i?/^ 



cos 



sin ^(/j^ / 



2aE 5i? as 



cos^?ag+ . 2 .a^fcdQ' 
+ sm 17- 



sm 



1? 9?? 



as d\ 



-sin^ ^?/2Q^ 



/2 



Pc 



+ sin2^/J>s(-T+ 



2 a2Ti+ 2cosi9aTi+ 2 / S^Tj 



OS 



siv? "d dip^ ~'~ sin ■!? 9?? ~'~ sin y d'ddip sin Sy? 



So 



So 



C Linearised matching in Poisson gauge for e = 

Here, as an example, we shall write the linearised matching conditions for the particular 
case of a flat e = FLRW region in the Poisson spacetime gauge in FLRW, for which 
W = x = Ya = 0. 

The odd part equations are given by 



2 ' 



W2-GI- —[U2 - (i?2+R^vr)lso] 

d 



So 



-aQay^ 



W2' - 2asGg - - 2aE(R5+i?;;,l'r)|So] = "^^^ 

Q2 - p2 - 2aE(i?5+R«,yr)] + 2/c"n^2 - {Rl+R'LVr 

The even part equations read 

T+ - [Wi - F'^ 

Qi+Yr- pi - 2aj:{R^+R^Y; 



3m 

.1 So 



a 



dg 

dr ' 



+ 2f~\Ui - iP^I^P^Lvn] + — = -/c'«s^, 

as or 



1 d 
as «A 



as (wi - - ^[Ui - (p,?+p;?.yr)lso]) 



So 



2m 



So 
AgV 



ciA aA V as 

So 2aE^,^,,(l) A 'D^ f25^ • .3^'^ 



as/c - 2m 
as "A V 



(58) 
(59) 

(60) 

(61) 

(62) 
(63) 

(64) 



+ 



2m 



as/c - 2m dr 



•9 . ^ 2as 
aj:—AeV + 



So 



as/c — 6m 5C/'^^^ 



3m 

as/c as/c - 2m 



(as/c - 2m)' 



r(2C/«- AgP), 



(65) 



(66) 



(67) 
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Jc as asjc 



1 



1 



2 So 



1 5Q ^ 2m d f5Q 



fc as a^fc V as / a|/| fc{ai:fc - 2m) 



1 dV as fc — 4m „ ^ 1 cos 6 d 



/c sin e de 



6Q_ ^ arji 



as as 
1 



/c 5r fc{aj:fc - 2m) 



2V + a^fc 



+ 2^(^ - ^p^yr) + yR- /c^ = fc { -2m— 



cos 



edv 



dV 



dr sm9 86 



dr J 



(68) 



(69) 



D Identification with Mukohyama's perturbation variables 

The variables N, fC and KL here (see subsection I3.2p have been chosen so that they correspond 
to those of [23] and [25], while the function r there corresponds to the Schwarzschild radius 
here, and takes the value as/c on the FLRW side of Sq. 

In the FLRW region, the S"^ scalar variables used here are related to the variables in 
[231 [25] in the following way: 



1=0 

oo 

Y^h^^Y = a\WR + W) 



1=0 

oo 



J2ht^Y = a\-2<^ + x) 

1=0 

oo ^ ^ 

J2 htY)Y = a2[-2$/2 + -//V + -As2X - -fx" + 2ff'YR + As2iUi + Ji) + H] 



1=1 

oo 



i>i 



1=1 

oo 



1>1 



h+^^Y = a\Q2 + 3^2' - 2^3^2) 



1=1 

00 



i>i 



Y h^L)Y = a\x' + Qi+Yr + yi' - 2^3^i - ^x) 
1=1 ■' J 

00 

y=a2(Z^2 + 3^2) 



1>1 



1=2 

00 



l>2 



E^(LL)^=«'(^l+3^1+oX) 



1=2 



l>2 



where Y stands for the corresponding spherical harmonic for any given pair (/, m), and where 
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the sum over the values of m is to be understood. Also note that, as explicitly indicated, the 
expressions on the right hand side must be restricted to the corresponding values of /. 

At either (it) side on Sq, the decomposition of the first fundamental form perturbation 
tensor q^^^ in \23\ I25j. is related to the Hodge scalars decomposition used here by 

oo 

1=0 

oo 

1=1 

oo 

1=1 

oo 

Y.'^iLT)Y= R\^2 (70) 

1=2 

oo 

1=0 

oo 

J2'^iLL)Y= P\;,2- 
1=2 

Analogously, the relations regarding the second fundamental form perturbation k^^^ follow 
from the above replacing o" by k on the left hand side, and q hy k on the superscripts on the 
right hand side quantities. 

Finally, the perturbation of the vector Z at either side is represented in [25] by 

oo oo oo oo 

1=0 1=0 1=1 1=1 

where Q stands for the Q used in [25] . 



References 

[1] Einstein A and Straus E G "The influence of the expansion of space on the gravitation 
fields surrounding the individual stars". Rev. Mod. Phys. 17, 120-124 (1945); "Correc- 
tions and additional remarks to our paper: The influence of the expansion of space 
on the gravitation fields surrounding the individual stars", erratum ihid 18, 148-149 
(1946) 

[2] Mars M and Senovilla J M M "Geometry of general hypersurfaces in spacetime: junction 
conditions". Class. Quantum Grav. 10, 1865-1897 (1993) 

[3] Krasiiiski A (1997) Inhomogeneus Cosmological Models, Cambridge University Press, 
Cambridge 

[4] Bonnor W B "A generalization of the Einstein-Straus vacuole" , Class. Quantum Grav. 
17, 2739-2748 (2000) 

[5] Fayos F, Jaen X, Llanta E and Senovilla J M M "Matching of the Vaidya and 
Robertson-Walker metric". Class. Quantum Grav. 8, 2057-2068 (1991) 



28 



[6] Fayos F, Senovilla J M M, Torres R "General matching of two spherically symmetric 
spacetimes", Phys. Rev. D 54, 4862-4872 (1996) 

[7] Oppenheimer J R and Snyder H "On continued gravitational contraction", Phys. Rev. 
56, 455-459 (1939) 

[8] Cocke W J "Some collapsing cylinders and their exterior vacuum metrics in General 
Relativity", J. Math. Phys. 7, 1171-1178 (1966) 

[9] Shaver E and Lake K "Nonstationary general-relativistic 'strings' ", Phys. Rev. D 40, 
3287-3301 (1989) 

[10] Senovilla J M M and Vera R "Impossibility of the cylindrically symmetric Einstein- 
Straus model", Phys. Rev. Lett. 78, 2284-2287 (1997) 

[11] Mars M "Axially symmetric Einstein-Straus models", Phys. Rev. D 57, 3389-3400 
(1998) 

[12] Mars M "On the uniqueness of the Einstein-Strauss model", Class. Quantum Grav. 
18, 3645-3663 (2001) 

[13] Nolan B C and Vera R "Axially symmetric equilibrium regions of Friedmann-Lemaitre- 
Robertson- Walker universes", Class. Quantum Crav. 22, 4031-4050 (2005) 

[14] Mena F C, Tavakol R and Vera R "Generalisation of Einstein-Straus models to 
anisotropic settings", Phys. Rev. D 66, 044004 (2002) 

[15] Hartle J B "Slowly rotating relativistic stars: I Equations of structure", Ap. J. 150, 
1005-1029 (1967) 

[16] Chamorro A "A Kerr cavity with a small rotation parameter embedded in Friedmann 
universes", Gen. Rel. Grav. 20, 1309-1323 (1988) 

[17] Gerlach U H and Sengupta U K "Junction conditions for odd-parity perturbations on 
most general spherically symmetric space-times", Phys. Rev. D 20, 3009-3014 (1979) 

[18] Gerlach U H and Sengupta U K "Even parity junction conditions for perturbations on 
most spherically symmetric space-times", J. Math. Phys. 20, 2540-2546 (1979) 

[19] Martfn-Garcia J M and Gundlach C "Gauge-invariant and coordinate-independent per- 
turbations of stellar collapse II: matching to the exterior", Phys. Rev. D 64, 024012 
(2001) 

[20] Cunningham C T, Price R H and Moncrief V "Radiation from collapsing relativistic 
stars. 1. linearized odd-parity radiation", Ap. J. 224, 643-667 (1978) 

[21] Cunningham C T, Price R H and Moncrief V "Radiation from collapsing relativistic 
stars. 2. Linearized even-parity radiation", Ap. J. 230, 870-892 (1979) 

[22] Carter B and Battye R A "Gravitational Perturbations of Relativistic Membranes and 
Strings", Phys. Lett. A B35, 29-35 (1995) 

[23] Mukohyama S "Perturbation of junction condition and doubly gauge-invariant vari- 
ables". Class. Quantum Grav. 17, 4777-4798 (2000) 



29 



[24] Mars M "First and second order perturbations of hyper surf aces" , Class. Quantum Grav. 
22, 3325-3347 (2005) 



[25] Mars M, Mena F C and Vera R "Linear perturbations of matched spacetimes: the gauge 
problem and background symmetries", Class. Quantum Crav. 24, 3673-3689 (2007) 

[26] Dolezel T, Bicak J and Deruelle N "Slowly rotating voids in cosmology". Class. Quan- 
tum Crav. 17, 2719-2737 (2000) 

[27] Czinner V and Vasuth M "Revisiting rotational perturbations and the microwave back- 
ground". Int. J. Mod. Phys. D 16, 1715-1723 (2007) 

[28] Vera R "Symmetry-preserving matchings" , Class. Quantum Crav. 19, 5249-5264 (2002) 

[29] Stephani H, Kramer D, MacCallum M A H, Hoenselaers C and Herlt E (2003) Exact 
solutions of Einstein's field equations. Second Edition, Cambridge University Press, 
Cambridge 

[30] Stewart J, "Perturbations of Friedmann-Robertson- Walker cosmological models". 
Class. Quantum Crav. 7, 1169-1180 (1990) 

[31] Noh H and Hwang J "Second order perturbations of the Friedmann world model", 
Phys. Rev. D 69, 104011 (2004) 



30 



